Methane Hydrate under High Pressure: 
Searching for Centering of Hydrogen Bonds 
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The structural, electronic, and spectroscopic properties of a high-pressure phase of methane hy- 
drate (MH-III) are studied by first principles electronic structure calculations. A detailed analysis 
of the atomic positions suggests that ionization of hydrogen-bonded water molecules occurs around 
40GPa and centering or symmetrization of hydrogen-bonds occurs around 70 GPa. These pressures 
are much lower compared with ioninzation around 55 GPa and centering around 100 GPa in pure 
ice. The transition may be observed with low-temperature IR/Raman spectroscopy of OH stretching 
modes or neutron diffraction. 

PACS numbers: 62.50.+p,82.75.-z,71.15.Pd 



Methane hydrate (MH), known as Burning Ice, is a 
special class of ice that contains methane molecules in 
cages or networks of hydrogen bonded water molecules. 
Low pressure phase of methane hydrate (MH-I) forms si 
structure of cages 0. MH-I, abundant in the deep ocean, 
has been attracting attention of the industry as a key ma- 
terial of new energy resource, whose amount is estimated 
twice as much as the total fossil fuel reservejlj. 

MH is also known as an important material for un- 
derstanding the mystery of the atmosphere of Titan, the 
largest satellite of Saturn. The conventional theory |2( 
could not explain abundant methane gas in Titan's at- 
mosphere because MH-I inside Titan was assumed to de- 
compose into ice and methane around 1 or 2 GPa, and 
escape to the atmosphere to be photodecomposed in early 
stage of Titan's history. To understand this mystery is 
one of the goals of the Cassini-Huygens spacecraft, which 
started its journey in 1997 and will arrive at Saturn sys- 
tem in 2004 Q. 

On the earth, in 2001, Loveday et al. (3, discov- 
ered new phases of MH by X-ray and neutron diffraction 
experiments under high pressure: MH-I transforms to 
MH-II ( sH cage structure) at 1 GPa, and then to MH-III 
phase (orthorhombic filled ice structure) at 2 GPa, which 
survives at least up to 10 GPa. Other researchers re- 
ported similar high pressure phases[H@]- Recently Hirai 
et al. |3 reported that MH-III survives up to 42 GPa at 
room temperature. Shimizu et al. Q have measured the 
site- and pressure-dependence of CH- and OH-vibration 
frequencies in these phases up to 5.2 GPa. Discovery of 
these high pressure phases allows us a new explanation of 
abundant methane gas in Titan's atmosphere: methane 
gas may be reserved in thick layers of MH-III under Ti- 
tan's surface and gradually emit to the atmosphere from 
the reservoir Q. 
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In this article, we focus on the features of MH- 
III as a promising material for investigateing center- 
ing or symmetrization of hydrogen bonds between wa- 
ter molecules. Advantage of using MH-III for study- 
ing the centering is that it is expected to occur at 
much lower pressure than that of p ure ice (ice VITice 

x transition) [13 EH Ep, E3, [13, Ea EH 113, EH dp, I2O1, 

which may make the difficult high-pressure experiments 
easier. In the following, we calculate the crystal structure 
and vibrational spectra of MH-III by using the density 
functional theory, so that we can predict and analyze ex- 
perimental results. 

We modeled the crystal structure of MH-III at 3 GPa 
by using the Pmcn symmetry, the lattice parameters 
and the position of atoms determined by the diffraction 
experiment p|. We chose the orientation of methane 
molecules carefully so that the molecules become close 
packed. The model consists of four methane and eight 
water molecules in the unit cell, of which one methane 
and two water molecules are symmetrically inequivalent. 
This structure was found stable after full geometrical op- 
timization, in which the enthalpy H = E + PV was min- 
imized by varying the lattice vectors and the positions 
of atoms without any constraints such as crystal symme- 
try. The structures at higher pressures ( Figs. ^ ) were 
calculated in a similar manner. 

The details of the electronic structure calculation |2lJ 
are as follows: the valence wave functions are expanded 
in a plane wave basis set truncated at a kinetic en- 
ergy of 1520 eV. The electron-ion interactions are de- 
scribed by the Vanderbilt-type ultrasoft pseudopotentials 
[22j. The effects of exchange-correlation interaction are 
treated within the generaliz ed g radient approximation of 
Perdew et al. (GGA-PBE) [23J. The Brillouin zones are 
sampled with 4x2x2 Monkhorst-Pack k-points [24| by 
using time-reversal symmetry only. In the geometrical 
optimization, the total stress tensor is reduced to the or- 
der of 0.01 GPa by using the finite basis-set corrections 

m 

Figure |3 shows the pressure dependence of the cal- 



FIG. 1: Crystal Structure of MH-III (a) 40 GPa, (b) 60 GPa, (c) 80 GPa 




FIG. 2: Lattice parameters: The solid lines indicate the 
present density functional theory calculation. The symbols 
indicate the experimental results of Hirai et al. || 



FIG. 3: Three symmetrically different distances between oxy- 
gen atoms, d(00), as a function pressure: the dashed and 
solid horizontal lines indicate the d(OO) at which ionization 
and centering occur in pure ice, respectively |T^|. 



culated lattice parameters at the zero temperature. It 
agrees well with the experimental data at room temper- 
ature by Hirai et al. Q. The cell is the most soft along 
the c-axis. This feature is also evident from the Figures 
1. The compression of the cell along the c-axis is mainly 
caused by the flattening of the two graphite- like wrinkled 
sheets normal to c-axis. As a result, the hydrogen bond 
network at very high pressures becomes sp2-bke structure 
in contrast to sp3 structure of ice Ih. 

In pure ice, centering of the hydrogen bond O — H - 
occurs as the oxygen- oxygen distance d(00) decreases 
with pressure increase [Hj, |n|. At large d(00) the 
proton occupies one of two potential minima along the 
O — H ■ ■ ■ O bond. As d(0O) becomes smaller the barrier 
between the two minima becomes lower. When d(00) 
becomes smaller than some critical distance, the proton 
occupies the single minimum at the midpoint between 
the two oxygen atoms (centering). In the following, we 
neglect the effect of finite temperature and quantum na- 
ture of hydrogen atom unless otherwise stated. Figure |3 
shows d(00) as a function of pressure. There are three 
curves corresponding to the three symmetrically inequiv- 
alent O - H-0 bonds. Around 70 GPa d(00)'s be- 
come sufficiently short for the hydrogen bond centering 
m pure ice Indeed Figures show that the center- 
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FIG. 4: Hydrogen position: The ratio d(OH)/d(00) calcu- 
lated by density functional theory is shown as a function of 
pressure. 



ing occurs around this pressure. Figure0|shows the ratio 
d(OH) / d(00) as a function of pressure. The ratio starts 
from 0.35 (water molecule) at 3 GPa and reaches to 0.5 
(centering) around 70 GPa. 

Now, let us look into the thermal effects. Benoit et al. 
|12| recently proposed three-stage scenario for explain- 



3 



— 3GPa 
■ ■ ■ 40GPa 

— 80GPa 




0.5 

d(OH)/d(00) 



FIG. 5: Distribution of hydrogen atom in a hydrogen bond of 
MH-III at 300K. 
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FIG. 6: Normal mode frequencies calculated with the density 
functional linear response theory as a function of pressure. 



ing the hydrogen bond centering [13| of pure ice with 
increase of pressure at room temperature: ice stays in 
molecular state under low pressures where H2O remain 
water molecule, then hydrogen atoms start to jump be- 
tween two potential minima in ionized state under mid- 
ium pressures, finally hydrogen atoms sit at the midpoint 
between the oxygen atoms in centering state under high 
pressures. The d(OO) at which ionization and centering 
occur in pure ice are indicated dashed and solid horizon- 
tal lines in Fig. |3| From this figure we can read that the 
ionization and cengering in MH-III are expected around 
40 GPa and 70 GPa, respectively. The probability distri- 
butions of d{OH)/d{00) for a hydrogen bond of MH-III 
at room temperature and with pressure 3 GPa, 40 GPa, 
and 80 GPa are calculated with Car-Parrinello molecu- 
lar dynamics 26] (Fig. |3J) , which clearly shows that three 
stage senario is also valid in MH-III but with much lower 
pressure than pure ice. 

Figure [fj] shows the pressure dependence of normal 
mode frequencies of MH-III calculated with the density 
functional linear response theory [27|]. The modes be- 
tween 3000 cm -1 and 3500 cm -1 are the CH vibrations, 
whose frequencies monotonically increase with pressure. 
The modes between 2000 cm -1 and 3000 cm -1 are OH- 
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FIG. 7: Vibrational spectra of MH-III at (a) 3 GPa and (b) 
80 GPa calculated by Car-Parrinello Molecular Dynamics are 
compared with the normal mode frequencies calculated with 
the density functional linear response theory. Solid line in- 
dicates result of CPMD at 300K, dashed line at 30K, and 
crosses are normal mode frequencies calculated with the den- 
sity functional linear response theory [2^|. 



stretching modes. Their frequencies decrease monoton- 
ically up to 70 GPa, where centering occurs, and then 
start to increase. This tendency is qualitatively in accor- 
dance with that of pure ice 14]. 

In Figure [7{a), vibrational spectra of MH-III at 3 GPa 
are calculated as a Fourier transform of velocity-velocity 
correlation function obtained from Car-Parrinello molec- 
ular dynamics j2|| and compared with the normal mode 
frequencies calculated with the density functional linear 
response theory 27].. The peak around 3100 cm -1 char- 
acteristic to OH-stretching vibration of hydrogen bonded 
water molecules is prominent at 30 K but disappears at 
300K. This fact is consistent with the experimental ob- 
servation that the Raman peak around 3100 cm -1 dis- 
appeared in MH-III at room temperature , and can be 
interpreted as the result of strongly anharmonic poten- 
tial. In Figure db) shows vibrational spectra of MH-III 
at 80 GPa. The OH-stretching modes of centered hy- 
drogen bonds are located between 2000 cm' 1 and 3000 
cm' 1 . At low temperature the peaks agrees well with the 
normal modes, while the peaks are blue-shifted at room 
temperature probably due to weak anharmonicity. The 
dissapearance of the peak of the OH-stretching mode in 
molecular and ionized states at room temperature and 
the survival of the peak due to OH-stretching mode in 
centered hydrogen bond in MH-III are analogous to those 
in ice VII and ice X [l7j . 

Figure |H1 shows the phonon dispersion of MH-III at 
80 GPa calculated with the density functional linear re- 
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FIG. 8: Phonon dispersion of MH-III at 80 GPa calculated 
with the density functional linear response theory |27f . 



(MH-III) with the density functional theory, and showed 
that ionization and centering of hydrogen bonds in MH- 
III may be observed around 40 GPa and 70 GPa, re- 
spectively, which are much lower pressure than those of 
pure water ice. Therefore MH-III may provide precious 
information about ionization and centering of hydrogen 
bonds between water molecules as the second example 
after pure ice. From the view point of planetary science, 
studying physical and chemical properties of clathrate 
hydrate such as MH will become more and more impor- 
tant to interpret and to understand the information of 
the outer solar systems sent from space-missions such as 
Cassini-Huygens spacecraft arriving at Saturn system in 
2004. 



sponse theory |27j. All frequencies at all wave numbers 
have positive frequencies indicating the mechanical sta- 
bility of this structure at zero temperature. However, 
further experimental and theoretical studies are neces- 
sary to prove the thermal stability at room temperature. 

In summary, we studied the structural and spectral 
properties of high pressure phase of methane hydrate 
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